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ABSTRACT 


The report outlines progress in implementing and refining two near-wall 
turbulence models in which the near-wall region is divided into either two or 
three zones. These models were successfully applied to the computation of 
recirculating flows. The research was further extended to obtaining 
experimental results of two different recirculating flow conditions in order 
to check the validity of the present models. Two different experimental 
apparatuses were set up; axisyrnmetric turbulent inpinging jets on a flat 
plate, and turbulent flows in a circular pipe with an abrupt pipe expansion. 

It is shown that generally better results are obtained by using the present 
near-wall models, and among the models the three-zone model is superior to the 
two-zone model. 
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NOMENCLATURE 



E 

H 

h 

k 

L 



Nu 


P 


P 



r 


T 

T 

w 


U 

U 


u-^ 


u, V, w 


V 

X 

y 

+ 

y 


coefficients in turbulence model 

nozzle diameter or diameter of pipe downstream from expansion 
empirical constant in logarithmic law 
nozzle-to-plate distance or step height (= (D - d)/2) 
heat transfer coefficient 
turbulent kinetic energy Cu^/E) 

nozzle length or length of the pipe upstream from expansion 
Nusselt number 

turbulence energy generation rate 
pressure 

Reynolds number based on diameter of pipe downstream from 

expansion/nozzls diameter 

radial coordinate 

temperature 

wall temperature, 

fluid temperature 

mean velocity in x direction 

friction velocity { =>/t / p) 

dimensionless velocity (U/U^) 

turbulent fluctuating velocities 

mean velocity in r direction 

coordinate parallel to flow 

distance from wal 1 

dimensionless distance (- yU /v) 

t 


2 


’“’eff 

e 

K 

M 




a 


effective diffusivity 

dissipation rate of turbulence energy (= v( 
von Kartnan constant 

dynamic viscosity 

effective viscosity (= y + 

turbulent dynamic viscosity 

turbulent kinematic viscosity (= y^/p) 

density 

Prandtl number 


2 


T 

i> 


turbulent Prandtl numbers for diffusion of k, e, 
and temperature 

turbulent shear stress 
dependent variable 


Subscripts 

B 

E, N, S, W 
e, n, s, w 


F 


k, e 


P 

t 

V 

w 


values at the edge of buffer layer 

values at east, north, south, and west node points 

values of the cell boundaries--east, north, south, 
and west sides 

values in the fully turbulent region 

values pertaining to kinetic energy and dissipation rates, 
respectively 

values at node point P 

turbulent values of quantity 

values at the edge of viscous sublayer 

wall values 
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1 . INTRODUCTION 


Many real flows of engineering interest contain regions of highly 
turbulent flow which exhibit separated, r^^attached, and recirculating flows 
even though the flow field is nominally steady overall. Aerodynamic engineers 
are interested in the influence of these flows on the heat and momentum 
characteristics on the wall boundary of airfoils, fluid machinery, blades, 
etc. It is frequently observed that the disturbance in the main stream 
usually has a significant effect on the wall boundary. This fact is usually 
seen in the situation of the flow in a pipe with a sudden expansion in which 
the separated flow reattached on the larger pipe wall and causes high 


heat/mass transfer rates. As another example, if we look at the turbulent 
jets when these impinge on a flat plate or body, it is also observed that the 
shear stress or heat transfer coefficient on the wall within the impingement 
region is significantly influenced by the turbulence intensity in the jet. 


This second example depicts the influence of the disturbance in the main 
stream on the heat, mass, and momentum characteristics on the wall. 

Therefore, considerable care must be taken for the evaluation of the 
wal 1-proximity region in handling the computational method. Generally, the 
near-wall region is divided into a viscous-affected .egion and a strongly 
turbulent region. Although the thickness of this near-wall region is usually 
two or more orders of magnitude less than the overall width of the flow, its 
effects extend over the whole flow field since, typically, 50% of the velocity 
change from the wall to the main stream occurs in this region. 

Generally, there are two methods of accounting for these wall-proximity 
region in numerical methods for computing turbulent flow: the wall-function 
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method and the low-Reynolds number modeling method. The former has some 
advantages: 1t is economical since computer time and storage required are 

much less than those required In the latter; it allows the Introduction of 
additional empirical Information; and it produces relatively accurate results 
by using fewer node points within the boundary layer compared with the 
low-Reynolds number modeling method, since the wall effect is evaluated only 
in the numerical cells next to the wall. However, the incorporation of 
complete wall-function in the wall adjacent cells requires enormous algebraic 
manipulation when many complicated equations are to be solved, such as the 
Algebraic Stress Model or Reynolds Stress Model. 

In this report the wall-function method was adopted to the k e 
turbulence model and near-wall models were developed for both the k and e 
equations. The reasons for the adoption of the k - e model are (1) it 
is generally recognized as a reliable model, (ii) it is relatively simple 
compared with multi-equation models of higher-order closures, and (iii) it is 
still used by many aerodynamic researchers since the k - e model 
produces results as good as those obtained using Algebraic Stress Models when 
applied to a compressible flow. 

Although the turbulence model is formulated carefully, sometimes 
predictions cannot be improved because of the numerical method used. The 
numerical model which is employed needs to be reviewed as well as turbulence 
models since the errors could be accumulated by numerical method per se. It 
is imperative to employ a numerical method which produces not only stable 
solutions but also accurate results. For example, the hybrid scheme of 
central and upwind developed by Gosman et al. [1] has been used by many 
researchers for a decade because of its excellent stability and simplicity 
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when applied to turbulent flows. Despite the advantages of this hybrid 
scheme, it cannot give accurate results especially near the region of cell 
tleynolds numbers being equal to 2. Near this region, the diffusive quantities 
are suppres.^sed which results in a large discrepancy between this predicted 
value compared to the analytical solution for a one-dimensional flow case. 
Consequently, the merit of higher-order closure turbulence models is easily 
shielded. Considering these points, a new numerical method was also developed 
in this report and was compared with the hybrid scheme. 

Many computations were made for two types of recirculation flows! 
turbulent impinging jets on a flat plate, and turbulent flows in a circular 
pipe with a sudden expansion. For the aim of validation tests of the present 
computational and turbulence model, experiments were performed for these two 
recirculating flows. 
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2, NUMERICAL METHOD 


2*1 Governing Equations 

The present work Is based on the numerical solution of the axisymmetric 
two-dimensional form of the time-averaged continuity, Navier-Stokes and the 
high Reynolds number version of k c turbulence equations. The 
equations following this approach for the present flow configuration can be 
written in the following general form; 

:pC*^(i^pU4) + -jjjr ( r pV (|>) 3 

= ( 1 ) 

where <j> stands for different dependent variables (U, V, k and s) for which 
the equations are to be solved. All the equations used in this report are 
summarized in Table 1. 

2 • 2 Numerical Solution Procedure 

The control volume approach was adopted for solving Eq. (1) in the finite 
difference scheme. The grid system used in this program is a so-called 
staggered grid system in which the value of each scalar quantity is associated 
with every grid node (i.e., the points where the grid lines intersect), 
although the vector quantities (velocity components) are displaced in space 
relative to the scalar quantities. This grid system has advantages in solving 
the velocity field since the pressure gradients are easy to evaluate and 
velocities are conveniently located for the calculation of convective fluxes. 

The finite difference scheme used in this paper is a method derived by 
expanding the exponential finite difference scheme of Spalding [2]. This 
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scheme can be ahown in the following form if Eq. ^1} is written in a finite 
difference form. 

Ap 4 >p =* b ( 2 ) 

where 

\ '^e'^ ^’•'e* 

Aw = '•^w'^ ■*■ ^''w* 

An =* Dp f{ iRpi) + max {-F^, 0} 

A$ = tig f( iRg i) + max {F^ , 0} 

Ap = Ag + Ajj + Aju + Ac 

b = S. fiVol 

9 

f( iR i) = max {0, 1 - iR i + -|^ iR iR i^} (3) 

D = X 

F = pU 

R = F/D = cell Reynolds number 

The function f ( iR i) shows the curve as in Fig. 1 in which the hybrid 
and the exponential schemes are compared with the present scheme. The nature 
of the hybrid scheme is such that it is identical with the central difference 
scheme for the cell Reynolds number range -2 5 R £2, and outside this 
range it reduces to the upwind difference scheme in which diffusion has been 
set equal to ^ero. However, as can be seen in Fig. 1, the departure of the 
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hybrid scheme from the exact solution (exponential scheme) is rather large at 
R '- + 2; also, it seems rather premature to set the diffusion effects equal to 
zero as soon as iR i exceeds 2- Considering the shortcomings of the 
exponential scheme, i.e., exponentials are expensive to compute and the scheme 
is not exact for two- or three-dimensional situations, the above expression of 
Eq. (2) is obtained by expanding the exponential expression to the fourth 
order term. This scheme is not particularly expensive to compute compared to 
the exponential scheme. Note that this scheme reduces to the upwind 
differencing for iR i greater than 4. Furthermore, accuracy has been 
improved in the range 1 < iRi < 4. 
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3, TURBULENCE NEAR-WALL MODEL 
3*1 Near-Wdll Two-Zone Model 

While viscous effects on the energy-containing turbulence motions are 
negligible throughout most of the flow, the no-slip condition at a solid 
interface always ensures that, in the immediate vicinity of a wall, viscous 
effects will be influential. Although the thickness of this viscous-affected 
zone is usually two or more orders of magnitude smaller than the overall width 
of the flow, its effects extend over the whole flow field since, typically, 

50 % of the velocity change from the wall to the free stream occurs in this 
region. 

The near-wall model which evaluates the mean generation rate and mean 
dissipation rate in the k-equation in the control volume cell adjacent to the 
wall was proposed by Chieng and Launder [3]. This model was applied to the 
computation of a turbulent impinging jet by Amano and Neusen [4]. In this 
model, as can be seen in Fig. 2, a parabolic variation of the turbulent 
kinetic energy is assumed which cornsonds to linear increase of fluctuating 
velocity with distance from the wall within the viscous sublayer. The 
turbulent kinetic energy, k, varies linearly towards the outer node points. 

t 

The turbulent shear stress is zero within the viscous sublayer, and the shear 
stress undergoes an abrupt increase at the edge of the sublayer while varying 
lintarly over the remainder’of the cell. The details of this treatment of 
k-equation are given in reference [3]. However, these local variations of 
turbulent quantities were not incorporated in the evaluation of both 
generation and destruction terms of the e-equation, but the value of c in 
the near-wall cell was approximated under local equilibrium conditions as: 
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where kp is the value of turbulent kinetic energy at the node point P 
adjacent to the wall. 

In the present study, the treatment of the c-equation in the near-wall 
cell is developed taking into consideration that the value of c near the 
wall is an order of magnitude larger than that in the fully turbulent core and 
reaches its maximum at the wall. Each term in the c-equation should be 
evaluated in accordance with the k-equation rather than being approximated 
under local equilibrium conditions. The procedure oP developing the two-zone, 
near-wall model is described herein. The two-zone model is defined to be the 
one in which the near-wall region is divided into two distinct regions: a 
viscous sublayer region (0 < y"** < 11), and a fully turbulent region 
(11 < y < 400). Now, noting that we have the relation near the wall 
C5]: 


or 


k « 3.5 

T 




U 

T 


If we define the Reynolds number based on the turbulence energy at the edge of 
viscous sublayer 




^ V 

li 


( 3 ) 


then Ry is found to be 20 which corresponds to y^ = 11.0 

In the viscous sublayer and in the fully turbulent region, turbulence 
energy, k, energy dissipation rate, e, and turbulent shear stress, t, are 
expressed as follows: 
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(i) Viscous sublayer 

k . k,(|-)' 


e 


2v (%~) 


(4) 


T = 0 


(ii) Fully turbulent region 


k = 



y + (kp - 



= by + a 


e 




(5) 


T = 




where 


and 


a 


= k - '‘P ' y 

'^P Vp - Vn ^P 



In the above equations, the notations in Fig. 2 are used. 

Ttie terms in the e-equation represent a e-balance involving 


convective transport, diffusion, generation, and destruction of c (see both 
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£q. (1) and Table 1). The convective and the diffusive terms are of minor 
influence near the wall and are neglected toward tiw wall* This fact is 
assured in the expression of k and e given in Eq* {4} as; 


which corresponds to no alffusion of k and e to the wall. The diffusional 
flux of energy dissipation out of the cell at its north, west, and east 
boundaries (in Fig. 2) are handled by the same differencing scheme that is 
employed over the remainder of the flow region. 

The mean generation and destruction rate in the e-equation can be 
obtained as follows: Over the fully turbulent region, the main velocity 
parallel to the wall is assumed to vary with distance from the wall according 
to 


Uk 


1/2 


Vp 



log (EVy^/v) 


( 6 ) 


* 1/4 * 1/4 

where k = kC andE = EC and k and E are the von Karman constants, 

U y 

The generation rate of k can be written as: 


P = + iV) 

' ay ax' 


where x = turbulent shear stress. 

Since the turbulent shear stress is zero within the viscous sublayer, by 
using Eqs. (5) and (7) the mean generation rate of e can be expressed as: 


(C 




(by a )^^^ 
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[t,., + ( 
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y 

) ](— + 
^ yn 


aU ^ aV 
lx 


) dy 


( 8 ) 
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By using Eq. (6) for the main velocity distribution, we can obtain the 
mean generation rate of e as: 


(C, f-p) 


, 1/2 


^r\ 




J/2 


(2(l<y^ - ky^) + ai)] + ^ [ 


n 


n 


T„i2(ky^ - '<y^) 




2 ^n " \ 1 3/2 

3 y^ b ^^n 


- k 


3/2 

V 



where 


/ 1 


ky-a'/2)(,y2.al«) 


X =< 


-1 


1/2 


JJ2 [tan"' (4) - tan"^ {-JL) ] 


1/2 


^ (-a) 


(a > 0) 


( 10 ) 


{a < 0) 


Unlike the generation rate of e, the destruction rate of e is not zero in 
the viscous sublayer. In the viscous sublayer, by using Eq. (4), the average 
value of k can be obtained as: 



and the mean destruction rate is expressed as: 
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Uv\ 
'^v •'n 


12 


2 ^n-^v 



( 12 ) 


vvhere Ry is defined in Eq. (3). 

To evaluate the destruction rate in the fully-turbulent region, we 
integrate (C«c^/k) between y and y . After including the linear 

« V n 

variation of k in Eq. (5) and the contribution of destruction in the viscous 
sublayer given in Eq. (12), the mean rate can be given as: 


"^2 k-) = 


12 

■i'n'^v 


^ 2 

(JL) 


1 - yv/^n 


2ab 




•^n ' 


log/+ b^) 


(13) 


where a and b are given in Eq. (5). The near-wall two-zone model described 
above is summarized in Table 2. 

3 . 2 Near-Wall Three-Zone Model 

In the previous section, we developed the approximation of near-wall 
region by dividing the region into two distinct zones. However, most of the 
experimental data show [6 and 7, for example] that both the linear and 
logarithmic profiles deviate from the experimental data in the buffer zone, as 
shown in Fig. 3. Therefore, in this section, a three-zone near-wall model is 

■f’ 

proposed which is comprised of a viscous sublayer (0 <y < 5) adjacent 
to the wall, a buffer zone (5 <y^ < 30), and a fully turbulent zone 
(30 < y"^ < 400). 

Figure 4(a) shows a computational node P whose associated control volume 
is bounded on the south side by a wall. In this figure the three zones are 
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shown such that the node point P lies outside the buffer zone assuming the 

near-wall cell is large enough. However, the computer program self adjusts 

the cases when P lies in either buffer zone or viscous sublayer. The behavior 

of the turbulent kinetic energy, k, and the turbulent shear stress, t, needs 

to be proposed especially in the viscous sublayer and the buffer zone. 

Bakewell and Lumley [8] reported that the streamwise fluctuating velocity 

increases linearly with distance from the wall, y, within the viscous sublayer 

and then it increases with y*^ in the outer region of the viscous sublayer. 

These data are plotted in Fig. 5. Hence, as we have done in the previous 

section, k is assumed to vary in a parabolic profile within the viscous 

suolayer. If we examine the slope of the fluctuating velocities in Fig. 5, 

the powers n of the fluctuating velocities in three directions fall in the 

range between 0.175 to 1.175 within the buffer zone. The average of the 

values n shows approximstely 0.7. However, m assumed that n is approximately 

^05 

0.5 which gives more insight by losing little accuracy, since u y * 
indicates k ~ y, i.e., linear approximation of k within the buffer zone. 

Within the fully turbulent region, the variation of k is controlled by the 
methodology of finite difference and the linear variation between node P and 
its northern neighbor is applied to extrapolate to the edge of the buffer 
layer (see Fig. 4(b) ) . 

Unlike the variation of k, the turbulent shear stress t hes a different 

dependence on y. If the streamwise velocity has a relation U - y near the 

wall, then, from the continuity equation, it is deduced that V ~y and 

2 

the convective acceleration is proportional to y . Hence, we have 
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This cubic profile generally yields very small values of -r within the 
viscous sublayer, which thus permits x to be treated as negligible in this 
viscous sublayer. While in the buffer zone, t is assumed to vary with cubic 
profile and undergoes a relatively sharp increase at the edge of the buffer 
layer and varies linearly over the remainder of the cell. This approximation, 
shown in Fig. 4(c), again gives more insight since, by setting t = 0 within 
the viscous sublayer, algebraic expressions can be simplified considerably 
with little loss of accuracy. 

From the above discussion, the variation of k, e and t can be 
summarized in the following form in the three zones. 


(i) Viscous sublayer 

k = k (J)‘ 
J'v 


e 


= 2v( 


3y 


1/2 2 


( 14 ) 


T = 0 


(ii) Buffer layer 
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(iii) Fully turbulent region 


k = 



^n -^B 


y + (kp - 


- h y 

yp - y^ 


P) 


= by + a 


c 


= k 



(16) 


where 



‘w 




yT^^P 


b = 


*^n - 
^n 


In Eqs. (14) - (16), the notations in Fig. 4 are used. 

The mean values of generation and destruction rates for both k and e 
equations are obtained by integrating the local variables of generation (P and 
C-jPe/k) and destruction (e and Cge /k) over the computational cell 
after inserting the relations in Eqs. (14) - (16) and then dividing by the 
volume of the cell. The streamwise velocity in the fully turbulent region can 
be given the following form. 


Uk 


w 


1/E 

B 

n 


_ 1 


-w 


In (E* y kg'^^/v) 


( 17 ) 


where E* and k* are defined in Eq. (16). 
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(i) The mean generation rate of the k-equation, P: 

The generation rate, P, can be approximated as; 

P = t(| + f) (18) 


Thus, 



By using Eq. (17) for the second term in Eq. (19), we can obtain the following 
form: 


■p = („4 . 4) (iU) 


172 — ^ " 7"^ 

PK*kB^ y„ 


ys« 






3V, 


( 20 ) 


in which (aV/sx)^ is assumed to be negligible. 
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(ii) The mean dissipation rate of the k-equation, T 

cp dy] (21) 

where the subscripts v, B, and F stand for viscous sublayer, buffer layer, and 
fully turbulent region, respectively. After substituting Eqs. (14) - (16) 
into Eq. (21), we obtain 




+ I (k„^^ - 1<B^^) + 2a(ky^ - kj^^) + a^x] 


where 


/ 1 

T72 


X 




(a > 0) 


1 <„ 1/2 


k. 1/2. 


^ (-a) 


Jan"' - tan"' (-^) ^ 


(a < 0) 


and where a and b are given in Eq. (16). 


(iii) The generation rate of the e-equation C^Pe/k: 

By using the expression in Eq. (18), we obtain as: 



dy 



£ 

k 




dy 


( 23 ) 



ORIGINAL PmU U 
OF POOR QUALITY 


20 


and substituting Eqs. (14) - (17), 


C ,lp « -iL r 2 ^ J /2 n jtj 


. 1/2 .. 1/2 


pkJ/2,^ 


J] 4. b .V 


(24) 


+ '\— ^ (2(ky^ - + ak)} 


^ (O.I.1/2 .,l/2i j. X 2 'n ■ ’'w 1 ,,.3/2 ,.3/2, , sVn 

+ ‘y 2 (k„ - kg )+ak)+-j.-jp — ,5 (k„ -kg )),jj] 


2/1, . 


(iv) The mean destruction rate of the e-equation, Cg t/k 


n 2 

e 


2 C„ / -^v 2 / -^B 2 / 

^2 f " ” 7 " M dy + I (|-) dy + I ( p ) dy] 

^ ^ sf /y^ B ^yg F 


(25) 


Since the k value has a singular point at y = 0, we approximate the 
integration within the viscous sublayer as we did in the two-zone model (see 
Eq. (11)). Then after substituting Eqs. (14) - (16), we obtain: 


2 19 kw 2 kr, 2 , y,. 

^2 F “ ^2 7T" 


( 26 ) 


^ ■ '^B'^'^n , a^ , 2ab 

+ 5 (--- - + 


‘Vn ^n " " ^B 


log Tp + b^)] 


All the results developed above are summarized in Table 3. 
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4. EXPERIMENTAL APPARATUS AND PROCEDURE 


4.1 Turbulent Jets Impinging on a Flat Plate 

Experimental heat transfer coefficients were obtained by using an 
electrlcally'-heated flat plate with an axisymmetric air jet Impinging normally 
to the surface of the plate (F1g. 6). Nozzle-to-plate spacingi distance from 
the nozzle centerline, and nozzle Reynolds number were varied so that the heat 
transfer characteristics of this configuration could be determined over a 
range of conditions. 

The test section was located in a large transparent enclosure to ensure 
that air movement within the room did not affect the tests. Holes in the top 
of this enclosure allowed heated air to escape; thus, the temperature within 
this enclosure did not vary appreciably during tests. The electrical power 
used to establish the uniform wall heat flux boundary condition was supplied 


by a DC power supply (1% ripple). The air flow rate was measured with a 
calibrated rotameter. Air and wall temperatures were measured with 
copper-constantan thermocouples. Thermocouple voltages, test-section voltage 
drop, and the voltage drop across the calibrated shunt, which was used to 
measure the current flow through the test section, were measured with a 
digital voltmeter. 

The test section consisted of a thin (0.81 mm) stainless steel plate with 
the heated portion being approximately 15 cm by 15 cm (see Fig. 7). Copper 
bus bars were soldered to both ends of the plate to ensure a uniform voltage 
drop along the test section. The electric resistance heating resulted in a 
uniform wall heat flux. The plate was mounted on a transite block, and its 
back and sides were heavily insulated to minimize conduction losses. Five 
thermocouples were mounted to the back of the plate 2.5 cm apart along one 
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diagonal of the test section and centered at the intersection of the two 
diagonals. A sixth thermocouple was located 5.0 cm from the center of the 
plate on the other diagonal. This arrangement permitted the symmetry of the 
jet to be checked by examination of the heat transfer coefficients. 

The nozzle consisted of a circular stainless steel tube 4.6 mm I.D, and 
L/0 = 75. A smooth plastic tube with a slightly larger I.D. than that of the 
nozzle connected the flow meter to the nozzle. The long nozzle length ensured 
a fully-developed velocity profile of the nozzle exit. The jet air 
temperature was measured with a thermocouple which was located in the plastic 
tube approximately 50 cm from the entrance to the nozzle. 

Before the electrical power was applied to the test section, the nozzle 
was positioned perpendicularly above the plate, at a specified H/D value at 
the intersection of the two diagonals where a thermocoupls was located. The 
symmetry of the resulting flow was checked in two ways. Prior to heating the 

plate, a velocity impact probe was used to measure jet velocities about 3 mm 

above the plate surface and 2,5 and 5 mm from the nozzle centerline along the 
two diagonals. This procedure confirmed a generally symmetric flow pattern, 
but indicated a slight instability in the flow which caused a small random 

shifting of the flow pattern. The symmetry was also checked by examining the 

heat transfer coefficients. This inspection also confirmed the symmetric flow 
pattern. After data were taken at one position, the nozzle was shifted 
laterally for a distance of one nozzle diameter along the diagonal and 
additional data were taken. 

Heat transfer data were taken after the power was turned on and steady 
state conditions (as indicated by temperature measurements) were obtained on 
the test section. For each test run, the air flow was set to a specified 
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level and was kept at this value through the run# The flow rate was very 
steady. Test section power was set to a level such that the temperature 
difference between the wall thermocouple at the jet stagnation point and the 
jet temperature was about 16K. The power was kept at this value throughout 
the runs at a set Reynolds number and H/D ratio. 

The heat flux from the test section wall was calculated using the voltage 
drop across the test section wall and the current flowing through it. lo 
calculate the outer wall temperature, utilizing the experimentally measured 
inner wall temperature and the heat generation rate, the steady-state 
differential conduction equation for a plane wall was solved, assuming 
one-dimensional heat flow and constant thermophysical properties of the 
steel. Radiation from the test section was estimated and found to have a 
negligible effect. Heat transfer coefficients were calculated u the wall 
heat flux, outer wall temperature and jet temperature. 

Data were taken at H/D ® 4, 7, and 10, Reynolds numbers at 10^, 

2 X 10^, and 4 x 10^ and -12 < r/D <22. Heat fluxes of 1950 to 8200 
w/m were used, which resulted in temperature differences (T^ - T^:) of 15K 
to 75K. Uncertainties in the experimental quantities are estimated to be; 
q", + R6 q, + 3%; and h, + 5%. 

4.2 Turbulent Flows in a Circular Pipe with an Abrupt Pipe Expansion 

To determine the heat transfer characteristics downstream of a circular 
abrupt expansion, three test sections were tested over a wide range of 
Reynolds numbers. The test sections were placed vertically in a flow loop, a 
schematic which is shown in Fig. 8. Liquid R-113 was pumped through the 
test sections at predetermined flow rates by adjusting the amount of fluid 
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allowed to bypass the flow loop. The inlet temperature to the test section 
was controlled by adjusting the electric preheater power. High pressure 
maintained single-phase flow throughout the tsst section. The flow was 
measured with calibrated rotameters. Fluid and wall temperatures were 
obtained with copper-constantan thermocouples. 

The test section (See Fig. 9) was constructed from a 12.2mm I.O., 457mm 
long stainless steel tube. A uniform wall heat flux was obtained by passing a 
DC current through the tube. The wall tnickness of 0.254mm helped to minimize 
axial conduction. Brass bus bars were silver soldered to each end of the 
tube, A thick layer of insulation reduced heat losses to a negligible level, 
Thi ‘ty thermocouples were spaced along the outside tube wall, including three 
on the inlet bus bar to help evaluate the effect of the axial conduction. The 
first 19 thermocouples were spaced 4.76mm apart; the next four were spaced 
12.7mm apart; and the final seven were spaced 25.4mm apart. The shorter 
thermocouple spacing just downstream of the expansion was necessary since 
relatively large variations over short distances in the heat transfer 
coefficients were anticipated in this area. The thermocouple voltages as well 
as the test section voltage drop and snunt voltage drop (from which test 
section current was obtained) were measured with a digital voltmeter with an 
accuracy of Ijiv. 

The nozzles were machined from a hard plastic ("Delrin"). Nozzle 
diameters of 2.38mm, 4.76mm, and 7.15mm were used which resulted in expansion 
ratios (d/0) of 0.195, 0.3D1, and 0.586, respectively. The nozzle exit was 
perpendicular and concentric to the stainless steel tube. Total nozzle length 
was about L/d 50 for all three cases so that a fully developed turbulent 
flow would be ensured at the nozzle exit. The entrance to the 
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nozzle was contoured so that there was a relatively smooth transition from the 
25.4mm supply piping to the nozzle. 

Before any data were taken, the flow rate and inlet fluid temperature 

were adjusted to give a specified Reynolds number based on the tube diameter 

downstream of the abrupt expansion. The test section power was then adjusted 

so that a maximum temperature difference between the wall and the fluid never 

exceeded approximately 35K, At the location of the maximum heat transfer 

coefficient, these power settings resulted in a minimum temperature difference 

4 • S 

ranging from about 2K at Re « 10 to lOK at Re « 1.5 x 10 . The 
temperature rise of the fluid ranged from about 3 to 5K. After steady-state 
conditions were attained, the data were taken. 

The inside wall temperatures were obtained using this measured outside 
wall temperature and the electrical power dissipation rate. The steady-state 
heat conduction equation with internal heat generation in cylindrical 
coordinates was solved by assuming one-dimensional radial conduction with 
constant thermophysical properties. Because of the thin tube wall, low 
thermal conductivity, and relatively high heat transfer coefficients, axial 
conduction was assumed to be small. Thus, the heat flux was calculated by 
dividing the total power dissipation by the inside wall area of the heated 
tube. The heat transfer coefficients were then calculated by dividing the 
heat flux by the local temperature difference (T^ - T^) at any location 
along the test section. 

Heat balances comparing electrical dissipation with enthalpy rise of the 
fluid generally were within ±S%. A propagation-of-error analysis suggests 
that the uncertainty in the heat transfer coefficients at the point of the 
maximum heat traitsfer coefficient range from about +4.7% for the largest d/0 
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ratio and highest Reynolds number to about +9.55J for the smallest d/D ratio 
and lowest Reynolds number. Uncertainties at other locations in the tube are 
smaller than these. The uncertainty in the Reynolds number is estimated to be 
± 2 %. 

5. PRESENTATION AND DISCUSSION OF RESULTS 

The experimental results of the heat transfer coefficients along the flat 
plate for the turbulent impinging jets are shown in Figs. 10-12 for three 
different nozzle Reynolds numbers, Roq, and for three nozzle-to-plate 
distances, H/D. These data show smoothly decreasing heat transfer 
coefficients with increasing distance from the stagnation point along the flat 
plate. The heat transfer coefficient also increases with increasing Re^ and 
decreasing H/D. 

As comparisons of these data with the numerical model, the near-wall 
two-zone model developed in Section 3.1 was used for the computation of the 
heat transfer coefficients along the flat plate for these impinging jet 
experiments and were compared with the experimental data, These computed 
results were further compared with the results obtained by employing the 
simpler near-wall models: the two-zone model of Chieng and Launder [3] in 
which the two-zone is considered only for the k-equations and not for the 
s-equations, and the one-zone model. These are shown in Figs, 13 and 14 for 
H/D « 4 and 10, respectively. For both cases the nozzle Reynolds number is 2 
X 10^. 

Note the generally close similarity of the results for the two different 
H/D; in particular, the calculated results show a much steeper slope in the 
impingement region. Comparing the one-zone model and the simpler version of 
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the two-zone model, the results are generally improved by approximately 30?^ by 
employing the two-zone model. Further comparison between the simpler two-zone 
model and the present two-zone model in which both the k and ^-equations are 
evaluated by taking the local values into consideration, it is obvious that 
the present effort improves the prediction, especially in the impingement 
regions approximately by 10/*5. However, the levels of heat transfer 
coefficients do not change in most of the regions except near the stagnation 
points. This indicates that the treatment of the e-equation is 
significantly sensitive, especially in the low Reynolds number region where 
very fine scale eddy motions take place. Consequently, such fine scale 
motions directly influence the behavior of e. 

Next the three-zone model developed in Section 3.2 is compared with 
experimental data of turbulent flows in a circular pipe with an abrupt 
expansion. The experimental results of this flow are shown in Figs. 15-19. 
Figures 15-17 represent the variation of Nusselt number for different Reynolds 
number based on a large pipe, Re^, and Figs. 18 and 19 represent the 
variation of Nusselt number for different pipe expansion ratios, d/D. For any 
particular test, the Nusselt number for different pipe expansion ratios, d/D. 
For any particular test, the Nusselt number starts at a low level, increases 
as one moves away from the pipe expansion section and then decreases, 
eventually reaching the fully developed condition in 30 to 40 pipe diameters. 
The maximum Nusselt number occurs about 6 to 8 step heights, H, downstream 
from the expansion section in all cases. The levels of Nusselt number 
increases with increasing Reynolds number and decreasing expansion ratio, d/0. 

Figures 20-23 display comparisons of the two-zone and three-zone models 
with the experimental data for two different Reynolds numbers and for d/D = 
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0.391 and 0,586. From these results, it is easily seen that the results 

obtained by employing the three-zone model show surprisingly better results 

than those obtained with the two-zone model. Note also that by using the 

two-zone model the predicted maximum Nusselt numbers display 5 - 205!^ higher 

values than the experimental data for ROq = 2 x 10^, while these 

predictions are 50 - 1005^ higher than the experimental values at Re^ = 8 x 

10^ when the two-zone model was used. On the other hand the predictions of 

the maximum Nusselt number with the three-zone model shows differences of -20% 
4 4 

at ReQ = 2 X 10 and +13% at Roq = 8 x 10 . Generally, it can be said 
that, by employing three-zone models, relatively better agreement with the 
experimental data can be obtained for a wide range of Reynolds numbers, but 
the two-zone model should not be used when the Reynolds numbers are greater 
than, say, 5 x 10^, since it predicts values for the heat transfer rate 
which are too high. 

Unlike the difference in the predicted levels of Nusselt number, the 

location of the predicted maximum Nusselt number does not vary with the 

different wall function models. Thus, it may be concluded that none of the 
various aspects of the wall function models developed in this study has any 
significant effect on the axial location of the maximum Nusselt number; the 
position of the maximum Nusselt number is controlled by the turbulence model 
per se. It is also concluded that the wall function treatment is not very 

C 

important for relatively low Reynolds number flows (Re^ < 5 x 10 ) but 
is significant for high Reynolds number flows. 

6. SUMMARIZING REMARKS 

The report has presented two different near-wall models developed based 
on the wall-functions. The principal advantage of these models, particularly 
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the three-zone model, is that the;/ produce wall proximity characteristics more 
accurately independent of the flow Reynolds numbers when applied to the 
recirculating flows. There are also several advantages noted in the use of 
these near-wall models developed here; 

1. a large number of cells is not necessary to obtain relatively 
accurate results if the near-wall model is employed; and 

2. we can improve predictions when applied to the k - e 
model with relatively little effort, ttc. 

It has been shown that although the three-zone model makes better 
predictions than the two-zone model, the latter still considerably improves 
the accuracy of prediction when compared with the one-zone model in which any 
local variations of turbulent quantities are not taken into account. 
Furthermore, modification of computer programs by incorporating the two-zone 
model is simpler than that required for the three-zone model. Therefore, when 
the computer program is very complicated, an adoption of the two-zone model 
can still improve the results. 

The studies of the application of these near-wall models to the flow 
around an aerodynamic wing by employing the k - e model and the Reynolds 
Stress Model will be presented in the next report. 
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Table 2. Near-Wall Two Zone Model - Cont. 
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Table 3. Near-Wall Three Zone Model 


Generation rate in k-equation P 


~ 4 . y 4)(JU) 


T (t - T ) 
n w' 

P«*Kb\ 


(1 




Dissipation rate in k-equation e 

I 

^',v y^c/3-’B 

+ I (k^^ - kg^) + 2a(k^^ - k^^) + a^x] 


Generation rate in e-equation (C^Pe/k) 


1 hilfi k i (1 _ (_y.)2} (M) 

7 C,y„ h ^yj ^ ^3y^B 


— — J— [t (— -.— + -X) 


+ J - kg*) + ax)] 






«f <*• 




« tM'kW*- • #•»♦ #«»■ ■• » 


ORIGINAL PAGE 
OF POOR QUALITY 


Table 3. Near-Wall Three Zone Model - Cont, 
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Fig. 3 Experimental verification of the inner-, outer-, 
and overlap-layer laws relating velocity profiles in 
turbulent wall flow 
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Experimental Set Up of Turbulent Impinging Oet on a Flat Plate 
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Fig. 10 Experimental Results of Heat 
Transfer Coefficient along a Flat Plate for 
Turbulent Impinging Jets 
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Fig. 12 Experimental Results of Heat 
Transfer Coefficient along a Flat Plate 
for Turbulent Impinging Jets 
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Fig. 13 Heat Transfer Coefficients along 
a Flat Plate for Turbulent Impinging Jets 
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Fig. 14 Heat Transfer Coefficients along 
a Flat Plate for Turbulent Impinging Jets 
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Fig. 18 Experimental Results of Nusselt Number Distribution for Turbulent 
Flows in Pipes with an Abrupt Pipe Expansion 
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Fig. 21 Nusselt Number Distribution for Turbulent Flow in a Pipe with an 
Abrupt Pipe Expansion 
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Fig. 22 Nusselt Number Distribution for Turbulent Flow in a Pipe with an 
Abrupt Pipe Expansion ‘ 
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Fig. 23 Nusselt Number Distribution for Turbulent Flow in a Pipe with an 
Abrupt Pipe Expansion 


